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ABSTRACT 


X-ray astronomy allows study of objects which may be associated with compact objects, i.e. neutron 
stars or black holes, and also may contain strong magnetic fields. Such objects are categorically 
non-spherical, and likely non-circular when projected on the sky. Polarization allows study of such 
geoemetric effects, and X-ray polarimetry is likely to become feasible for a significant number of 
sources in the future. A class of potential targets for future X-ray polarization observations is the 
high mass X-ray binaries (HMXBs), which consist of a compact object in orbit with an early type 
star. In this paper ws show that X-ray polarization from HMXBs has a distinct signature which 
depends on the source inclination and orbital phase. The presence of the X-ray source displaced from 
the star creates linear polarization even if the primary wind is spherically symmetric whenever the 
system is viewed away from conjunction. Direct X-rays dilute this polarization whenever the X-ray 
source is not eclipsed; at mid-eclipse the net polarization is expected to be small or zero if the wind 
is circularly symmetric around the line of centers. Resonance line scattering increases the scattering 
fraction, often by large factors, over the energy band spanned by resonance lines. Real winds are not 
expected to be spherically symmetric, or circularly symmetric around the line of centers, owing to the 
combined effects of the compact object gravity and ionization on the wind hydrodynamics. A sample 
calculation shows that this creates polarization fractions ranging up to tens of percent at mid-eclipse. 


1. INTRODUCTION 


High mass X-ray binaries (HMXBs) are among the 
brightest X-ray sources in the sky. They consist of an 
accreting compact object (a neutron star or black hole) 
in orbit with an early type star (O or B star). They 
represent an important stage in the evolution of binary 
stars with early-type components, and they dominate the 
X-ray output of galaxies with young stellar populations. 
HMXBs emit X-rays via accretion onto the compact ob- 
ject by gas from the strong wind from the companion 
star, or from Roche lobe overflow. The X-rays interact 
with the accretion flow and stellar wind, producing ob- 
servable signatures in timing and spectra. This provides 
a means to study the hot star wind and the accretion 


primarily at hard X-ray energies owing to obscuration, 
and transient sources which are detected primarily via 
long time-baseline monitoring. The most thoroughly 
studied systems include the objects Cyg X-1, Vela X- 
1, and Cen X-3. The majority of HMXBs contain pul- 
sating neutron stars. If so, orbital eclipses, X-ray pulse 
timing and primary radial velocities 
provides strong constraints on the orbital separa- 
tions and masses. Simple estimates suggest that the ac- 
cretion rate onto the compact object can be provided by 
the wind from the companion star in < half the sources, 
while Roche lobe overflow is required in the remainder 
(Contil197[Whital1985¢ [Kaper 1998). Though they are 


outnumbered by low mass X-ray binaries in our galaxy, 
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flow and the interaction of X-rays with these structures. 
One consequence of this interaction is the production of 
polarized X-rays via scattering. Polarization provides 
unique information into the geometry of the stars, wind 
and scattering region. In this paper, we explore the po- 
larization signatures associated with this interaction. 
There are ~10 HMXB sources with properties of their 
stellar winds and binary orbit which are well enough 
determined for quantitative modeling of the interac- 
tion of the accretion flow with X-rays (Kaped [1998). 
The number of known HMXBs is many times greater 
when sources with less well constrained properties are 
included ; these in- 
clude sources discovered by INTEGRAL which are seen 


these sources dominate the X-ray output from galaxies 
with larger star formation rates (Fabbiand|2006). Their 
mass transfer can be very rapid; in a system in which 
the primary has a radiative envelope and fills its Roche 
lobe, the primary radius will shrink as a consequence 
of mass transfer, but if the mass transfer is conservative 
the Roche lobe will shrink more rapidly, potentially lead- 
ing to a common envelope phase. If the primary is an 
evolved star with a convective envelope this outcome is 
more probable (Taam & Sandquist[2000). Common en- 
velope evolution when the primary is an evolved star is 
likely to lead to ejection of much of the envelope of the 
primary, leaving a binary consisting of a neutron star 
and the core of the primary. Whether due to nuclear 
evolution or to unstable mass transfer, their evolution- 
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ary lifetimes must be short (< 10° yrs). 

The wind from the companion absorbs and scatters the 
X-rays from the accreting compact object in HMXBs, 
and the column densities traversed by the X-rays range 
from ~ 1071 — 1074cm~?. The X-ray source luminosities 
range from ~ 10°° — 10°° erg s~+. In sources where the 
X-ray source is luminous and the wind is weak, the wind 
is ionized throughout most of the region between the two 
stars. In systems with more massive winds and weaker 
X-ray sources X-ray-ionization effects are a perturbation 
on the wind properties. Key questions about HMXBs 
include the nature of the compact object: its mass, vari- 
ability and intrinsic radiation pattern. Also of interest 
are the properties of the companion star wind and its uni- 
formity in density and temperature. The X-ray source 
can provide a probe for the study of wind regions which 
are only weakly affected by ionization or the gravity of 
the compact object. 

Many of the X-ray properties of HMXBs are affected 
by the stellar wind and accretion flow. This includes the 
variability around the orbit, which shows gradual eclipse 
transitions due to photoelectric absorption in the wind 
(Clark et_al] 1988), enhanced absorption at late orbital 
phases due to a wake or stream trailing the compact ob- 
ject (2006). 
The wind and accretion flow also provide torque to the 
compact object and regulate its spin or angular momen- 
tum. The structure of the wind and accretion flow is 
uncertain; winds from comparable single O or B stars 
have radiatively driven winds (Lamers & Cassinelli 
which are affected by instabilities (Owocki et_alJ[1988) 
shocks (Lucy & White/1980) and clumping (Cohen et al] 
[2006) ; observations provide the most reliable means for 
understanding these processes, via fitting of parameter- 
ized models. In addition, the compact object affects the 
gas flow, via its gravity and also X-ray heating and ion- 
ization. 

The next new astrophysical window will be the advent 
of measurements of X-ray polarization 
. Among other things, po- 
larization allows for potentially sensitive tests of the ge- 
ometry of astrophysical sources, on scales which are far 
too small to be imaged directly. There is only one source 
in the sky whose X-ray polarization is known, the Crab 
nebula (Weisskopf ot alll1978). Owing to visibility con- 
straints, it is likely that the first astronomical X-ray 
polarimetry observations will be of objects which have 
never before been observed with this technique. This 
motivates thorough and accurate modeling of the polar- 
ization properties of the brightest and (otherwise) best 
understood X-ray sources, for use as calibrators and test 
sources for X-ray polarimetry. HMXBs are among the 
sources best suited for this. Their orbital elements are 
relatively well understood and their orbital variability 
provides a predictably changing view with respect to an 
important source of polarization: the strong stellar wind 
from the companion star. X-ray polarimetry provides a 
means to probe the structure and the physical processes 
occuring in and near compact objects, on length scales 
too small to be imaged directly. In spite of scanty as- 
trophysical detections so far it is 
of interest to consider the possible signals and diagnos- 
tic use of X-ray polarization in known classes of cosmic 
sources, as it is likely that instruments with improved 











sensitivity will become available eventually. 
Polarization properties of HMXBs have been discussed 


previously in the context of scattering of optical and 
[Brown ot allligra 1978). These authors 
predicted that the linear polarization has a characteris- 
tic variability around the binary orbit, associated with 
scattering in the stellar envelope which is distorted by 
the gravity of the compact object. This variability can 
be used as a diagnostic of the shape of the star, and also 
of the inclination of the binary orbit. These effects are 
observed at the predicted level in some systems; their 
absence in other cases suggests that the optical light 
may be affected by contributions from structures whose 
shape is not symmetric around the line of centers. Re- 
lated work has been carried out b (1999); 
(2009): (2014). Useful 
formalism for the calculation of polarization has been 
provided by (1996). The polarization prop- 
erties of the accretion columns in X-ray pulsars have been 
calculated by (Meszaros et_al[1988). Predictions of X- 
ray polarization by the large scale structures in binaries 
has been carried out for distorted stars (Angell [1969); 
for Compton scattering within the accretion column of 
magnetic cataclysmic variables (McNamara et_al./2008), 
following up on pioneering work by (2004). 

X-rays have an important advantage over optical light 
which is that the source is almost certainly compact com- 
pared with the size of the binary system. That is, we ex- 
pect X-rays to be emitted from a region smaller than the 
Alfven shell, defined as the distance from the compact 
object where the pressure due to a dipole magnetic field 
balances the ram pressure of accretion. For fields ~ 10! 
Gauss this region has a size ~ 10° cm (Lamb et_al./1973). 
Far outside this region the influence of the field on the 
dynamics is negligible. The Alfven shell size can be com- 
pared with the primary star radius and orbital separa- 
tion which are ~ 10/7 cm. This disparity in sizes greatly 
simplifies the geometrical considerations associated with 
the use of X-ray polarization since the X-ray source is 
effectively a point source when considering the effect of 
the scattering by the stellar wind. Calculating the po- 
larization from these structures is straightforward, given 
certain simplifying assumptions, though it is important 
to take into account the effects of atomic absorption and 
resonance line scattering, which in turn depends on the 
gas dynamics, along with electron scattering. 

The polarization properties of HMXBs will be affected 
by the intrinsic polarization of the compact X-ray source, 
and by the polarization imprinted by the stellar wind 
and accretion flow. These can be distinguished by their 
differing variability behavior: the compact source gen- 
erally varies on a pulsation timescale for sources con- 
taining a pulsar, which is ~ seconds — minutes. The 
very few sources which contain black holes also have 
their strongest intrinsic variability on comparably short 
timescales. Wind variability, on the other hand, is as- 
sociated with the orbital timescale (~ days) or possibly 
on the wind flow timescale which is >hours. In some 
HMXBs the X-ray source is luminous enough to ionize 
the wind almost completely, so the light observed during 
and near eclipse, and its polarization, is dominated by 
electron scattering. 

In this paper we present a general discussion of the be- 


havior of X-ray polarization from HMXBs as a function 
of the parameters of the system and the viewing posi- 
tion. We present calculations of the polarization signa- 
ture for various simple analytic models for the wind den- 
sity and ionization structure, and discuss the dependence 
of these on parameters: wind optical depth, inclination, 
and viewing direction or orbital phase. We also explore 
the effects of wind hydrodynamics, i.e. departures from 
spherical symmetry due to the effects of the X-ray source 
gravity and heating. We utilize sample three-dimensional 
dynamical models. We take the intrinsic polarization of 
the compact object as unpolarized and explore the com- 
bined effects of geometry and scattering physics on the 
predicted linear polarization in the 0.1 - 10 keV X-ray 
band. 

In section we present numerical calculations for 
spherical winds using only electron scattering. The ef- 
fects of resonance line scattering are discussed in section 
[3] In section [5] we present models which include photo- 
electric absorption and an ensemble of resonance scat- 
terers. In section [6] we present models utilizing three- 
dimensional hydrodynamic models for the wind density 
and velocity field. We use these to derive approximate 
predictions for the polarization levels expected for sev- 
eral well known systems in section [7] The Appendix 
presents simple analytic estimates of the polarization for 
idealized HMXB conditions; many of these mirror earlier 
results. The computational techniques and level of de- 
tail we employ are similar to those of our previous work 
(Dorodnitsyn & Kallman||2010) on the polarization from 


warm absorbers in Syefert galaxies. 


2. SPHERICAL WINDS 


The basic geometry of an HMXB is illustrated in figure 
This shows a schematic of an HMXB with a typical 
orbital separation and a sample X-ray ionized zone. It 
also shows the geometry we use when calculating the 
polarization: the X-ray source is at the origin and the 
primary orbits in a plane inclined by an angle 7 to the 
line of sight. The orbital phase is described by the angle 
Oy relative to the line of centers. The system is viewed 
along the y axis. 

In this section we will consider the simple case in which 
the wind density is spherically symmetric about the pri- 
mary star, and the only interaction between the X-rays 
and the wind is electron scattering. If so, the scattering 
phase function follows the Rayleigh form 
i960), The scattered emissivity at each point depends 
only on the local gas density and on the flux from the X- 
ray source. We assume the X-ray source radiates isotrop- 
ically with a constant total luminosity Lo and is unpo- 
larized. 

We calculate the polarization signatures in the form of 
the three Stokes parameters for linearly polarized light 


using the formal solution to the equation of transfer 
(Mihalas|1978). 


L(e) 1 + cos? y 
{90} = [wee r)S(e,r)e~7©") ¢ sin? y cos(27) 
U(e) sin? x sin(27) 


(1) 
where ¢ is the photon energy, S(¢,r), is the source func- 
tion, K(e,r) is the opacity, r(e,r) = f K(e,r)d¢ is the 
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optical depth from a point r to a distant viewer, y is 
the scattering angle and 7¥ is the angle between the scat- 
tering plane and a reference direction on the sky. Here 
and in what follows we describe the Stokes parameters 
in terms of the luminosity seen by a distant observer, 
ie. the total energy (in ergs s-' sr~*) radiated by the 
system in that direction. Equation [I] defines the scat- 
tered luminosity, L (polarized plus unpolarized); obser- 
vations are also affected by an unscattered component, 
Ly(e) = Lo(e)e~7©"*) where rx is the position of the X- 
ray source. For the purpose of evaluating these quantities 
we adopt a cylindrical coordinate system with the X-ray 
source at the origin and the viewing direction along the 
y axis. The angle 7 corresponds to the azimuthal angle 
on the plane of the sky and is measured relative to a 
line which is the intersection between the orbital plane 
and the plane of the sky. In the case of pure electron 
scattering the source function is S = Lo/(4rr?) where 
r, is the distance from the X-ray source and the opacity 
is K(e,r) = n(r)orn. 

Some useful results are presented in the Appendix. 
The fractional polarization relevant to observation is 
P= /(Q? +U?)/(L+ L,). This quantity is zero when 
the system is viewed at inclination 7 = 7/2 and at either 
conjunction, i.e. orbital phase 0.5 or 1, Oy = 0 or a. 
The polarization at i = 7/2 is amaximum at quadrature, 
phase 0.25 or 0.75. We can also define the polarization of 
the scattered radiation only, i.e. P; = \/(Q? + U?)/L, 
and the value of this quantity at quadrature and i = 77/2 
depends only on the orbital separation and on the dis- 
tribution of gas density in the wind. The value of P 
is less than the value of P, by a factor proportional to 
the Thomson depth whenever the X-ray source is not 
eclipsed. At i = 0, i.e. when the system is viewed face- 
on, the polarization rotates with orbital phase at a rate 
twice the orbital rotation rate. 

As an illustration of these results, we calculate the 
polarization produced by a single-scattering calculation 
by numerically evaluating equation The wind veloc- 
ity is assumed to be radial relative to the star with a 
speed given by u(r) = vo + (Ucoo — Vo)(1 — rx /r) where 
the terminal velocity is vo>>=1000 km s~! and vg=100 


km s~!, and r, = 10! cm. The wind mass loss rate 
is Mwing = 107-?Moyr7!. The orbital separation is 
a= 1.5r. 


In the remainder of this paper we discuss quantitative 
calculations of polarization produced by wind scatter- 
ing. These include calculations based on both the very 
simple spherical wind with parameters given in the pre- 
vious paragraph and also on hydrodynamic calculations 
of the wind density structure which make no assump- 
tions about spherical symmetry. In the former case, it 
is illuminating to estimate the accretion rate onto the 
compact object if the accretion is supplied solely by the 


wind. This can be done using simple estimates based on 
Bondi & Hoyle (1944!) and [Davidson & Ostriker (1973), 


ie. Mace = TR?..NxVexMy Where the accretion radius is 


Ry = 2GM/vz, M is the mass of the compact object and 
Uz is the speed of the wind at the X-ray source. This cor- 


responds to Mace/Mwina =1.8.~x 10-4(M/Mo)*u5¢ a5 
where vzg = vz/(1000kms~*), a12 = a/(10'2cm) or an 
X-ray luminosity L, ~ 1.8 x 10% ergs! using the pa- 
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rameters given in the previous paragraph and assuming 
an efficiency of converting accreted mass into energy of 
0.1. This luminosity is less than the time averaged lumi- 
nosities of most HMXBs by factors ~5 — 100, which likely 
reflects the fact that wind mass loss rates may be greater, 
additional mass can be supplied by Roche lobe overflow, 
and the wind and accretion flow dynamics can be en- 
hanced by the influence of the compact object gravity 
and ionization. The results in what follows, those which 
are based on the simple spherical wind approximation, 
must interpreted subject to this caveat. 

An additional caveat which applies to essentially all 
of the results presented in this paper is the assumption 
of single scattering. The validity of this assumption de- 
pends on the wind optical depth from the X-ray source; 
multiple scattering is important when this quantity ap- 
proaches or exceeds unity. For parameters similar to 
those discussed so far, this quantity for a distant ob- 
server corresponds to Trp ~ Nea ~ 2 x 10-4 Mgajy vz 
where Mg = M/(10-8Moyr7'). Comparing this expres- 
sion with the accretion luminosity estimate given above 
implies that multiple scattering can be important when 
the accretion luminosity exceeds ~ 10°%ergs~+ for a1 Mo 
compact object. The likely effect of multiple scattering is 
to produce smaller net polarization than for single scat- 
tering, since the second and subsequent scatterings have 
a wider range of scattering angles and planes than for 
single scattering. If the depth is moderate, i.e. < 10, 
there will still be a significant fraction of photons which 
reach the observer after a single scattering. If so the net 
polarization will likely be less than predicted here, by 
factors of order unity, and our results should be modified 
to include multiple scattering effects for such sources. 

Maps of polarization projected against the sky are 
shown in figure[2] This shows contours of constant inten- 
sity as solid colors separated by solid black curves, along 
with lines corresponding to polarization vectors, which 
appear as dashed curves. These are plotted vs position 
in units of 10! cm. These illustrate many of the results 
presented in the Appendix: At high inclination, i = 7/2, 
and at orbital phase angles 0 and 7 (orbital phases 0 
and 0.5) polarization vectors are perfectly circumferen- 
tial, contours of constant intensity are also circular and 
the net polarization is zero. At orbital phase angles 7/2 
and 37/2 (orbital phases 0.25 and 0.75) the net polariza- 
tion is maximum. At inclination i = 0 the star always 
influences the polarization, so that the net polarization 
fraction is constant, but the position rotates with the 
orbit. 

More insight comes from the Stokes quantities as a 
function of the orbital phase. These are shown in figure 
[3]for inclinations i = 0 and i = 7/2. U is green while Q is 
shown in red; solid is 4 = 7/2 and dashed is i = 0. This 
further illustrates the results from the previous figure: 
at 7 = 0 the polarization oscillates between Q and U 
along orbital phase, and the vector sum is constant. At 
i = 7/2 U is always small (in this convention) and Q 
is a maximum at phase angles 7/2 and 37/2. Another 
way of displaying the same thing is shown in figure 
which shows U plotted vs Q, for 1 = O in green, and 
i = 7/2 in red. This shows that the trajectory of U vs. 
Q is circular at i = 0 and becomes linear at i = 7/2. As 


shown by (1978), this trajectory is an ellipse 


Kallman & Bautistal/2001 


for intermediate inclinations, and the cases in figureBlare 
the extremes of eccentricity. This has been suggested as 
a means for measuring inclination (Rudy & Kempli97a). 
Figure [5] shows the scattered polarization fraction vs. 
orbital phase, ie. P; = \/(Q? + U?)/L. This shows that 
the maximum polarization fraction is ~10% for the pa- 
rameters chosen here, when the polarized component is 
compared with the total scattered component. This is 
comparable to the result in the Appendix, calculated for 
a thin shell at the star rather than for an extended wind. 
We emphasize that this value depends primarily on geo- 
metric quantities: the extent of the wind, and orbital sep- 
aration relative to the stellar radius. It does not depend 
on the wind density or column density since it is a com- 
parison of scattered quantities. The difference between 
i = 0 and i = 7/2 is clearly apparent: the former pro- 
duces approximately constant polarization fraction and 
the latter oscillates between 0 and a value comparable 
to the i = 7/2 value. Figure [6] shows the polarization 
fraction measured relative to the total radiation, scat- 
tered plus direct, ic. ie. P = \/(Q?2+U7?)/(L+4+ Ly). 
This illustrates the diluting effect of the direct radiation 
for this model, which has tr, ~ 0.25. Maximum linear 
polarization at i = 7/2occurs just following the eclipse 
transition, where the departures from circular symmetry 
on the sky are not negligible, and where the direct radi- 
ation is blocked by the primary star. At mid-eclipse the 
polarization at i = 7/2 is zero due to symmetry. Fig- 
ure [7| shows the polarization angles from the same set 
of models. The angle sweeps through 180 degrees twice 
per orbital period for i = 0, while the angle is constant 
(though undefined near conjunctions) for 1 = 7/2. 


3. RESONANCE LINE SCATTERING: SINGLE LINE 


The models discussed in section [2] include solely elec- 
tron scattering; they do not include resonant scattering 
in bound-bound transitions. This process (when associ- 
ated with UV and optical transitions) is the dominant 
driving mechanism for the winds from early type stars. 
Here we consider scattering in the X-ray band. In the ab- 
sence of X-ray ionization, the ions which are most abun- 
dant in early-type star winds are in charge states from 
~1— 5 times ionized, and do not have many strong X- 
ray resonance line transitions. Ionization by a compact 
X-ray source can produce ions of arbitrary charge state, 
depending on the X-ray flux and gas density, and so res- 
onance scattering can affect the scattered X-ray inten- 
sity and polarization. The ionization structure in HMXB 
winds has been explored by 
who showed that the surfaces of constant ionization in a 
spherical wind are either nested spheres surrounding the 
X-ray source, or else open surfaces enclosing the primary 
star. The ionization distribution depends on a param- 
eter g = €/€, where &, = Lz/(nza*) where Lz is the 
ionizing X-ray luminosity, n, is the gas density at the 
X-ray source and a is the orbital separation. This is a 
particular example of the scaling of ionization in opti- 
cally thin photoionized gas, which depends on the ion- 
ization parameter £ = L,/(n(r)r2) (Tarter et_al][1969; 

. The properties of resonance 
scattering as applied to X-ray polarization in HMXBs 
are dominated by the fact that, owing to the restricted 
regions where ionization is suitable, most lines can only 
have strong opacity over a fraction of the wind. The 





Figure 1. Schematic of high mass X-ray binary (HMXB) as viewed from above the orbital plane, showing typical orbital separation and 
companion star size. Coordinates are labeled according to their use in the Appendix. 
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Figure 2. Map of intensity and polarization vectors projected on the sky for four orbital phases and two inclinations. Colors correspond 
to scattered intensity according to the color bar, in which the labeled quantities correspond to logig of the specific intensity. Short black 
lines correspond to polarization vectors. Spherically symmetric wind is assumed. 
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Figure 3. Stokes parameters for spherical wind, Q (red ) and U (green) vs. phase for inclinations 1 = 0 (dashed) and i = 7/2 (solid). 
The total intensity I is shown as dashed black for 1 = 0. At 1 = 0 the polarization oscillates between Q and U along orbital phase, and the 
vector sum is constant. At i = 7/2 U is always small (in this convention) and Q is a maximum at phase angles 7/2 and 37/2. 
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Figure 4. Stokes parameters for spherical wind, Q vs U i = 0 (green) and i = 7/2 (red). This again illustrates the circular path of the 
net polarization direction with orbital phase at i = 0 and the constancy of the direction at i = 7/2. 


regions where this occurs most often resemble spherical 
shells surrounding the X-ray source. 

Resonance scattering differs from electron scattering 
in its sensitivity to the velocity structure of the wind. 
At a given observed photon energy, scattering can occur 
only over a resonant surface with shape determined by 
the wind velocity law; in a spherical wind with monotonic 
velocity law these surfaces are open surfaces of revolution 
symmetric about the line of sight to the X-ray source. 

Resonance scattering affects polarization by redis- 
tributing the polarization among the components of the 
Stokes vector according to the phase matrix which is a 
linear combination of the Rayleigh and isotropic phase 
matrices; coefficients of the matrix depend on the angu- 
lar momentum quantum numbers of the initial and final 
states of the transition, 7 (Lee et_al/1994). The matrices 
describing this process have been calculated by{Hamiltonl 
(1947). In practice most of the strongest X-ray reso- 
nance lines have initial and final values jiower = 1/2 and 
Jupper = 3/2, which correspond to the Rayleigh phase 
matrix, and we adopt this in our calculations. Similar 
assumptions were employed in 
(2010). 

Polarization also depends on the scattering angle, and 
the resonance condition constrains the scattering geom- 
etry, thereby imposing a relation between the energy of 
the scattered photon and its polarization. Resonance 
scattering can have a cross section which is much greater 
than electron scattering, but only over the relatively nar- 
row energy band spanned by the resonance line, includ- 
ing the effects of Doppler broadening associated with the 
wind motion. Within this band, the polarization can be 
greater than that produced by electron scattering, owing 
to the greater cross section, given suitable geometry. 

In the remainder of this section we illustrate these ef- 
fects by generalizing the single scattering spherical wind 
calculations to include resonance scattering. We consider 
a single resonance line, chosen to crudely resemble a line 
such as O VIII La. The parent ion is assumed to exist 
over a range of ionization parameter 1 < log(£) < 3. We 
calculate the optical depth using the Sobolev expression 


(Castoi[1970): 


Ar 
v(r)(1+ x ea - 1) 
(2) 


where f is the oscillator strength, n(r) is the gas num- 
ber density, x;; is the ion fraction, y; is the element 
abundance, z is the position along the line of sight, and 
A is the line wavelength. For the source function we 
adopt the same expression used for electron scattering 
S = Lo/(4rr2) since thermalization is unimportant and 
the size of the X-ray source is small compared with the 
other length scales of interest for HMXBs. 

Results of our simple single line calculation are shown 
in figure [8] This shows the line profile as a function of 
orbital phase for a system viewed edge on (i = 7/2). For 
each phase we display the luminosity, with the trans- 
mitted luminosity shown in black and the scattered lu- 
minosity shown in green (bottom panel) and polarization 
fraction (top panel). The polarization angle is not shown 
because there is no significant dependence on energy or 


2 
Te 
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orbital phase when the system is viewed at high incli- 
nation; the polarization is always perpendicular to the 
orbital plane. The horizontal axis is energy in units of 
the wind terminal velocity. 

At phase 0 (superior conjunction of the X-ray source) 
the shape of the profile in luminosity is similar to P-Cygni 
profiles familiar from UV resonance lines in hot stars: the 
outflow produces blue-shifted absorption (shown as neg- 
ative energy in these figures) of the continuum. This 
absorption is offset from the zero of energy owing to the 
fact that the continuum source in this case is the compact 
X-ray source, and the wind speed at the X-ray source 
is approximately half the terminal speed. In addition, 
the X-ray source creates an ionized region within which 
line scattering cannot occur. As a result, the absorption 
occurs in a relatively narrow region of energy near the 
energy corresponding to the wind terminal velocity. The 
scattered emission, on the other hand, is essentially sym- 
metric in energy, since scattered emission comes from a 
region which is not necessarily along the line of sight to 
the X-ray source. Projection effects make the scattered 
emission appear at all energies |e| < equa./c where €o 
is the line rest energy. The scattered emission is unpo- 
larized at phase 0 for the same reason as in the pure 
electron scattering case. At phase 0.25 (quadrature) the 
wind is viewed perpendicular to its velocity vector at the 
X-ray source, and therefore the absorption covers a large 
fraction of the energy spanned by the wind. The scat- 
tered emission is polarized up to ~50%; the degree of 
polarization is symmetric around the center of the line, 
owing to the fact that the wind velocity structure is sym- 
metric around the line of centers. At phase 0.5 (inferior 
conjunction of the X-ray source) there is no transmitted 
flux (due to occultation by the star) and the emission is 
predominantly red-shifted owing to the location of the X- 
ray ionized zone in the receding part of the wind. In this 
case, even though the flux is all scattered, the polariza- 
tion fraction is negligible owing to the circular symmetry 
of the scattering region of the wind as viewed in the plane 
of the sky. 


4. ABSORPTION EFFECTS 


The ionization balance in an HMXB wind is deter- 
mined at each point by the X-ray flux and by the gas 
density. The X-ray flux in turn depends on the effects of 
geometrical dilution and on attenuation. In order to do 
so, we utilize results from XSTAR 
at each point in the wind in order to determine 
the opacity. This is done along radial rays originating 
from the X-ray source, so that the opacity at each point 
is calculated using an ionization parameter which corre- 
sponds to the flux transmitted along the corresponding 
ray. 

This single stream transfer treatment is similar to the 
transfer treatment used by XSTAR with one difference: 
XSTAR uses the local spectral energy distribution (SED) 
to calculation the ionization along a ray; here we use 
the same unattenuated SED everywhere for the ioniza- 
tion calculation (so we can use a stored table), but we 
calculate the ionization parameter self-consistently, i.e. 
by calculating the transmitted flux at each position and 
then using that value to calculate the ionization param- 
eter. This simplification allows for more efficient compu- 
tation without unduly sacrificing physical realism. This 
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Figure 5. Polarization fraction of scattered radiation i.e. Ps = \/(Q? + U?)/L for spherical wind, i = 0 (red) and i = 7/2 (black). 





0.1 


0.08 
T 
1 


fraction 














0 1 2 3 4 5 6 


phase angle 


Figure 6. Polarization fraction relative to total radiation, scattered plus direct, for spherical wind, i = 0 (red) and i = 7/2 (black). 


angle 





0 1 2 3 4 5 6 


phase angle 


Figure 7. Polarization angle for spherical wind, i = 0 (red) and i = 7/2 (black). 
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Figure 8. Spectrum, polarization fraction and angle in a resonance line vs. energy in units of the terminal velocity at orbital phase 0.25 
for inclination i = 7/2 
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approximation will be least accurate when the column 
density is highest, i.e. > 1073 cm~?. In such regions, the 
ionization is expected to be low, owing to attenuation, 
and X-ray scattering will be relatively unimportant. Our 
approximation will tend to produce less attenuation at 
a given column density than a self-consistent calculation 
would. 


5. RESONANCE LINE SCATTERING: ENSEMBLE OF LINES 


X-ray ionization creates an ensemble of resonance lines 
in an HMXB wind from the many trace elements such as 
C, N, O, Ne, Mg, Si, S, and Fe. At any point in the wind 
the ionization depends most sensitively on the ioniza- 
tion parameter, defined in the previous section. Figure[9] 
shows the distribution of ionization parameter produced 
in a wind with a density distribution which is spherically 
symmetric around the primary star, as discussed so far 
in this section. We adopt a spherical wind with mass loss 
rate 10-°Moyr—', corresponding to a total wind column 
density of 3 x 1073 cm~?, and an X-ray source luminos- 
ity of 10°° erg s~+. Each spatial region with a distinct 
ionization parameter will have a different distribution of 
ion abundances, and hence a unique distribution of res- 
onance line opacity. These can be calculated under the 
assumption of ionization equilibrium. Figure [10] show 
examples of the spectrum and the polarization fraction 
produced by such a distribution from the spherical wind 
shown in figure [9] Ionization and line opacities were 
calculated as a function of ionization parameter using 
the xsTAR code (Kallman & Bautistal2001)). These were 
binned into 1000 energy bins over energy from 0.1 eV 
to 104 eV. Since the Doppler shifts associated with the 
wind speed is less than the energy bin size, we adopt an 
approximate form for the source function each resonance 
line: 


Lo Veo — Vo 
Anr2 Aec/e 


where Tiine is the line optical depth from equation[2] and 
the quantity (Vo. — Up)/(Aec/e) takes into account the 
fact that the line does not conver the entire energy bin. 
This treatment assumes that the line optical depths are 
not large, which is an adequate approximation for our 
situation. Depolarization effects at large Sobolev optical 
depths associated with multiple scatterings are not taken 
into account. These calculations serve to illustrate the 
magnitude of the polarization effects expected from res- 
onance scattering in HMXBs. Simulations suitable for 
quantitatively diagnosing the wind or the X-ray source 
properties will require reexamination of these assump- 
tions. 

Results of spectra and polarization fractions are shown 
in figure[10] These demonstrate that the resonance lines 
cover a significant fraction of the X-ray energy band and 
that they can scatter a significant flux of X-rays, and 
create polarization fractions as high as 0.5, much greater 
than would be produced by electron scattering alone. 
A difference between resonance scattering and electron 
scattering is that resonance scattering in a given line oc- 
curs within a relatively narrow spatial region where the 
parent ion is most abundant. For X-ray lines, these are 
most likely to be approximately spherical surface sur- 
rounding the X-ray source. This region tends to pro- 


S(e,r) = (3) 


duce a small polarization which is almost independent of 
the viewing angle. Thus resonance scattering produces 
weaker modulation of the polarization with orbital phase 
than does electron scattering. 

It is also worth noting that the effects of scattering can 
be either polarizing or not polarizing, depending on the 
scattering angle. Also, lines which appear in absorption 
in the spectrum can have enhanced polarization in their 
troughs owing to the presence of scattered light in the 
residual intensity. In the spectrum shown in figure [10] 
the polarization in the absorption lines is generally lower 
than in the emission features, or in the adjacent contin- 
uum. This is due to the fact that the scattered light in 
the residual intensity is forward scattered, and so is not 
polarized because of the scattering angle. 


6. HYDRODYNAMIC MODELS 


The stellar wind in HMXBs is not spherically symmet- 
ric. Physical processes which affect the wind and ac- 
cretion flow in HMXBs include: the three-dimensional 
geometry of the flow (ie. the flow both in and out 
of the orbital plane), rotational forces, the influence 
of gravity and radiation pressure from both stars in 
the binary, transport of the radiation from the stars 
into the flow and the reprocessed radiation out of the 
flow, X-ray heating and ionization, and departures from 
thermal equilibrium due to advection and adiabatic 
heating and cooling. The dynamics of X-ray heated 
winds have been discussed b (1980), 
and [Hatchett & McCray (1977), and multi-dimensional 
models have been calculated b (1990); 
(1994); (2008 

In this section we illustrate the effects of the wind hy- 
drodynamics on the polarization. We use two sample 
wind hydrodynamic models calculated using a numeri- 
cal model similar to that described in 
but extended to three dimensions. The hydro- 
dynamic models are computed on one hemisphere of a 
spherical grid, assuming reflection symmetry about the 
orbital plane. A non-uniform grid of 448 (r) by 128 (0) 
by 512(¢) zones is used with highest resolution near the 
surface of the primary star and in the vicinity of the 
accreting neutron star. This model calculates the wind 
dynamics in three dimensions taking into account the 
radiative driving by the UV/optical radiation from the 
primary, and also the gravity of the star and the com- 
pact X-ray source. A fixed X-ray luminosity of 10°° erg/s 
is used to calculate X-ray heating using the approximate 
formulae given in Blondin (1994) and the effects of X-ray 
ionization on the dynamics via changes in the UV radia- 
tion force multiplier. The two models have identical di- 
mensions: primary radius 2.4 x10!? cm and orbital sepa- 
ration 3.6 x10!2 cm. They differ in their mass loss rates, 
which are 4 x10~’ Moyr7! and 1.7 x10-°Moyr—. The 
maximum wind speed in both models is approximately 
1600 km s~'. A plot showing the density contours and 
velocity vectors is shown in figure[L1] This clearly shows 
the influence of the gravity of the compact object in cre- 
ating a region of higher density and non-radial wind flow 
in the vicinity of the X-ray source. 

The distorted stellar wind of the hydrodynamic models 
is illustrated with velocity vectors in the orbital plane 
and a series of transparent density isosurfaces. Right 
panel is the model with a higher mass loss rate (Mdot 
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Figure 9. Distribution of ionization parameter in the orbital plane for a spherical wind with Muwind = 10-5 Moyr~?* illuminated by a 
10°8 erg s~! X-ray source. Contours are labeled with logé. 
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Figure 10. Spectrum and polarization fraction vs energy for spherical wind in units log(E/keV) at orbital phase 0.25 for inclination 
i=n7/2 
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= 1.7e-6); left panel is a lower mass loss rate (Mdot = 
4.0e-7). The lowest density isosurface corresponds to a 
density of 4e9 cm~? in both panels. Velocity vectors are 
not shown for values less than 800 km/s. 

The high mass loss rate model has a denser spherical 
component of the wind, but a less pronounced accretion 
wake. The low mass loss rate model has a larger volume 
of wind moving at low velocity due to photoionization of 
the wind in the vicinity of the X-ray source. The result 
is a denser wind coming off the primary along the line of 
centers of the binary system and a larger, denser accre- 
tion wake that wraps more tightly around the primary 
star. 

We apply the same calculation of X-ray scattering to 
the hydrodynamic wind as was done in section [3B] We 
assume an X-ray source luminosity of 10° erg s~! and 
a [=2 power law ionizing spectrum. The distribution 
of ionization parameter in the orbital plane is shown in 
figure [12] This shows a spiral structure, owing to a cor- 
responding structure in the gas density. 

Figures [13] shows the spectra and polarization pro- 
duced by the hydrodynamic wind model with M = 
4.0 x 10-’Meyr7! at orbital phases 0, 0.25 and 0.5. 
Comparison with figure[0]shows that the hydrodynamic 
model produces a spectrum which is similar to that from 
a smooth spherical wind. However, the degree of po- 
larization is significantly greater, due to the departure 
from spherical symmetry created by the X-ray source. 
In particular, the spectrum during eclipse is greater, ow- 
ing to the fact that the wake structure extends beyond 
the disk of the primary star during eclipse. The spher- 
ical wind produces fractional polarization averaged over 
the 0.1 10 keV energy band, of at most 10% during 
eclipse. The spherical wind produces polarization val- 
ues which are much less than this at mid-eclipse. The 
hydrodynamic wind produces fractional polarization of 
approximately 21% at mid-eclipse. 

Polarization calculations for hydrodynamic models 
have been carried out for the two mass loss rates shown 
in figure[I1] and for several choices of ionizing X-ray lu- 
minosity. Note that the hydrodynamic models assume a 
fixed X-ray luminosity (used to calculate X-ray heating 
and ionization within the simulation) of 10°° erg/s, in- 
dependent of the X-ray luminosity used to calculate the 
spectrum and polarization. Moreover, the prescribed X- 
ray luminosity is not generally consistent with the mass 
accretion rate derived from the hydrodynamic simula- 
tion, and thus these models are not fully self-consistent. 
Nonetheless, they serve to illustrate the behavior of the 
polarization and its dependence on wind density and X- 
ray luminosity. 

Values for polarization fractions during eclipse are 
shown in table [I] These are averages of the linear po- 
larization over energy in the range from 0.1 — 10 keV. 
Polarization values reflect competing effects. For highly 
ionized gas, the wind is essentially transparent and the 
polarization is due to Thomson scattering. The polariza- 
tion fraction is proportional to the Thomson depth (when 
small) and also to departures from circular symmetry of 
the scattering material in the plane of the sky. Partially 
ionized gas can produce larger polarization, owing to the 
greater cross section associated with resonance scatter- 
ing, although each line has very limited spectral range 


and the ensemble of lines for any given model seldom 
has a width which exceeds Ac/e ~ 0.1. On the other 
hand, partially ionized and near-neutral gas also is af- 
fected by photoelectric absorption. This limits the size of 
the X-ray scattering region, and tends to reduce the net 
polarization. Table[I]shows that, for the limited range of 
parameters spanned by our models, more highly ionized 
models tend to produce greater polarization. That is, the 
effects of photoelectric absorption offset the effects of res- 
onance scattering in partially ionized gas, leading to very 
weak dependence of the polarization on X-ray luminosity 
when the luminosity is not large. We expect that at very 
low X-ray luminosities this effect will be even stronger, 
since the partially ionized zone containing the resonance 
scattering gas will shrink and become more round, and 
photoelectric absorption will remove a larger fraction of 
photons. 


7. DISCUSSION 


Our results so far on the polarization produced by wind 
scattering in HMXBs can be summarized as follows: Po- 
larization depends on inclination: high inclination pro- 
duces variable polarization fraction plus constant angle; 
low inclination produces constant polarization fraction 
plus variable polarization position angle. The maximum 
attainable continuum polarization fraction scales approx- 
imately proportional to electron scattering optical depth, 
for depths less than unity; this is partly a consequence 
of our single scattering assumption and must be suitably 
modified if multiple scattering is important. At high in- 
clination, polarization fraction of the scattered radiation 
(Ps) at phase 0 and 0.5 is less than at phases 0.25 or 
0.75; for a spherical wind the phase 0 and 0.5 polariza- 
tion is zero. Polarization fraction of the total radiation, 
including unscattered, is small for the parameters consid- 
ered here for all orbital phases out of eclipse. A spherical 
wind thus produces narrow intervals of high polarization 
during eclipse but away from phase 0. Resonance line 
optical depths are greater than for electron scattering, 
and so can produce greater linear polarization. On the 
other hand, the optical depths are often large, and the 
X-ray scattering regions tend to be small, i.e. nearly 
circular around the compact object, producing smaller 
orbital phase modulation. The hydrodynamics of the 
interaction between the wind and the compact object 
breaks the spherical symmetry and increases the net po- 
larization. 

Predictions of the polarization signatures for partic- 
ular known HMXBs require hydrodynamic simulations 
for each system incoporating known physical parame- 
ters: the sizes and masses of the components, the primary 
wind mass loss rate and the X-ray source luminosity and 
spectrum. These would then need to be analyzed using 
a transfer calculation such as the single scattering cal- 
culations we have presented. In this paper we have pre- 
sented generic models for the wind densities and X-ray 
source properties, both for spherical winds and incopo- 
rating the hydrodynamic effects of the two gravitating 
centers. These generic models can be used to infer, in 
a very simple way, the polarization signatures expected 
from various sources. 

Table [2] shows the known parameters of the 5 bright- 
est and best-studied HMXBs. The relevant quantities 
are the source X-ray luminosity, primary type, mass loss 
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Figure 11. Density contour and vector plot of wind hydrodynamic models. Right panel is model with M=1.7x 10-6 Moyr~* left panel 
is model with M = 4.0 x 10-7 Meyr—!. Spatial scale is the same for both panels. Yellow contour corresponds to a density of 6 x 108cm~. 
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Figure 12. Distribution of ionization parameter in the orbital plane for hydrodynamic wind illuminated by a 1 erg s~! X-ray source. 
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Figure 13. Spectrum and polarization fraction vs energy at phase 0, 0.25, 0.5 for hydrodynamic wind with M = 4.0 x 10-7 Moyr! 
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Figure 14. Diagram showing estimated polarization fractions in lines and continuum at conjunction vs. quadrature for the objects in 


table 
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taken from [Conti (1978) and from{Kaper (1998). We do 


not include several systems which generally have smaller 
observed X-ray fluxes or less well constrained properties. 


Table 1 
Hydrodynamic Model Results. Values for 
the polarization fraction are given for 


several hydrodynamic models described With these quantities we can use our spherical wind 
in the text. models to crudely predict the polarization fractions and 
a orbital phase modulation of the polarization for these 
Model M Le Pg=0.5 sources. We do this in the following way: we use the 
Moyr-! erg s-! results in table[I]to construct an approximate scaling of 

1 4 x10-7 1038 0.25 mid-eclipse net polarization with the wind mass loss rate 

2 4 x10-7 10°” 0.25 and X-ray luminosity. We then apply this to the known 

3 4 x10~7 iy? 0.29 HMXBs in table 1. The results are shown graphically 

4 17 x10~8 10°8 0.22 in figure The most relevant quantities for each of 


the table 1 HMXBs, X-ray luminosity and wind mass- 

loss rate, are plotted on the axes and the values for each 

object are plotted as points. Contours of constant pre- 

Table 2 dicted mid-eclpise polarization are shown as solid curves, 

Paes Sens and F roper ics Neaded ist and labeled. This shows that high polarizations are ex- 

redicting Wind Polarization. : ; : 

pected for some objects, those with the strongest winds 

and weakest X-rays generally. This demonstrates that 


pblse ep; Be Mwina be polarization fractions in the range from 5 — 30 % are ex- 

Moyr71 erg s7! pected at mid-eclipse for these systems. More detailed 

Vela X-1 BO.5 Tab 7 x10~6 0.05 information is contained in the spectra and the time vari- 

Cen X-3° O6.5TI-IT 2 x1077 0.3 ation of the net polarization during eclipse; interpreta- 

ea oe ce a tion of these signals requires modeling which is tailored 
cael B06 lab 5 x10-7 3 to each particular system. 
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APPENDIX 
APPENDIX 


Illustrative analytic results can be derived for an idealized HMXB system: a spherical primary star with radius 
R ~ 10" em in a circular orbit with an X-ray source with separation a ~ 1.5R. We view the system from some 
phase angle Oy relative to the line of centers where Oy = 0 corresponds to superior conjunction of the X-ray source. 
In addition, the system can have inclination i defined such that 1 = 27/2 corresponds to viewing in the orbital plane. 
The X-ray source has a luminosity Lx ~ 10°” erg s~! and illuminates the primary star and wind isotropically. The 
wind has a mass loss rate M ~ 10-8Mo yr7!, a terminal velocity v.. ~ 108cm s~!, and a velocity law u(r) which 
depends only on the distance from the center of the star. If so, the effects of scattering in the wind depend on the 
wind column density, which can be characterized by the column density N = he n(r)dr and an approximate value for 


this is N = nza where ng = Milo is the density at the X-ray source, yz is the mean particle weight, and my is the 


hydrogen mass. For the fiducial parameters given above, nz = 3.2 x 108 cm~? Mgug ‘aj? and N = 3.2 x 107° cm~? 
Msgvg ‘aja. The ratio of the total intensity of scattered radiation to the total emitted radiation can be characterized 
by the quantity Lscate/Lx ~ No ~ 2 x 10-40 /ornMgvg ‘azo where ¢o is a characteristic scattering cross section and 
orn is the Thomson cross section. 

We consider a coordinate system in which the observer is along the y axis. Then the scattering plane is perpendicular 
to the (x,z) plane. We can rewrite equation [Jas: 


L n(n) 1+ cos? y 
Q?p= Loo f av = 4 sin” x cos(27) (1) 
U "t | sin? y sin(27) 





For the purposes of computation, a more convenient expression is in terms of cylindrical coordinates centered on the 
X-ray source, with the observer on the axis at infinity: 


L 27 oo co (p? + a0) 
= oO AE 2 cos 
pnteh aL mf acres [em 


where (p,y,¢) are the radial, angular and axial cylindrical coordinates with the axis along y. These are indicated on 
figure [I] In the case a star located at a distance a from the X-ray source, orbiting in a plane with a normal which is 
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inclined at an agle 7 with respect to line of sight, and with an orbital phase described by an angle 0,, this expression 
can be rewritten in terms of coordinates centered on the star: 


p=r((xcos 6, + (y +a) sin Oy)? + (zsini + (y + a) cos, — zsin@,) cosi)*)!/? (3) 
y= tan7! zsini + (y+ a)cos@, — aera) cosz (4) 
x cos 6, + (y+a)sin 6, 


¢ =—-r((y + a) cos 6, — «sin 6,) sini + zcosi (5) 


where x,y,z are cartesian coordinates centered on the star r is the distance from the center of the star and a = a/r. 
This can be evaluated analytically for special cases. 
In the case of a single star, a=0, 6,=0, and if the density distribution is axially symmetric,|Brown & McLean (1977 


have shown that 


L oo 1 2(1 + cos? @) + sin? i(1 — 3 cos? 6) 
Qe= Loo | ar | sin OdOn(r, 8) sin? i(1 — 3 cos? ) (6) 
U 0 1 0 


so that the fraction varies from zero for zero inclination to a maximum value determined by the departure of the 
density distribution from spherical. 

From equation [2] it is apparent that, in the case where the density is independent of 7, the net polarization is zero 
and the surfaces of constant polarization in the integrand are circular. This is expected to be the case for a binary 
system where the gas density is circularly symmetric around the line of centers, and when viewed at conjunction from 
the orbital plane (i.e. 6,=0 or 7/2 and i=7/2). In the case of a binary system viewed along the angular momentum 
axis, i.e. 2=0, the intensities can be shown to be: 


L Zo 
Q > = <2, cos 20, + Zo sin26, (7) 


U T, sin 20, + Zz cos26, 
where 

Lp = Loo i r-dr if aene\<—<—$—<— — — _. (8) 

0 0 (r? + 2ar sin @ sin? ¢ + a?) 

m DBD aps se ay ie Died 2 

T, = Lye | dr | ee ee (9) 

0 0 (r? + 2ar sin @ sin* ¢ + a?)? 

oo c oo Ogi 

Th= Loo | rar | sin aed = ee (10) 

) 0 (r2 + 2arsin 6 sin* 6 + a?)? 








and Zz = 0 by symmetry if the density distribution is independent of ¢. This clearly represents a circle in the (Q,U) 
plane. This, plus more general cases, were explored py nec all (1978). 

In the case of a binary viewed from within the orbital plane, it is straightforward to show that the polarization is zero 
at phases 0 and 0.5, i.e. superior and inferior conjunction of the X-ray source. At phase 0.25 or 0.75, i.e. quadrature, 
the stokes quantities can be written: 


4 





LY Loom [* (181? +150? — 5) patie — $F + (Zr — 6m ($28) 
Qp =e [ rdrn(r) § (2272 + 1902 + 2) calltars + SB — (2832 + 26)n (22) (11) 
0 


This can be evaluated analytically in the limit that the density distribution is a thin shell at r =r, in which case the 
net polarization Q/L varies from unity for a >> r, to 0.886 for a =r, 
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